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Abstract 

In this paper we present an estimate of accuracy for a piecewise 
polynomial approximation of a classical numerical solution to a non 
linear differential problem. We suppose the numerical solution u is 
computed using a grid with a small linear step and interval time T u , 
while the polynomial approximation v is an interpolation of the values 
of a numerical solution on a less fine grid and interval time T v «T U . 
The estimate shows that the interpolant solution v can be, under 
suitable hypotheses, a good approximation and in general its compu- 
tational cost is much lower of the cost of the fine numerical solution. 
We present two possible applications to linear case and periodic case. 



1 The problem 

Let [a, b], a < b, an interval on the real line R, and let C 1 ([0, +oo) x [a, b]) the 
space of continuously differentiate real functions denned on [0, +oo) x [a, b]. 
Let F : C 1 ([0, +oo) x [a, b]) — * R a continuously differentiable real functional. 
Then, if u G C 1 ([0,+oo) x [a, b}), the partial differential equation, usually 
associated to hyperbolic conservation laws, 



1 




du d , 

is in general non linear on the unknown u. Let % G C 1 ([a, 6]) and u a , two 
real numbers. We consider the boundary problem 

?(x) = 0Vi6(a,&) 

= u a , u(t,b) = u b Wt G [0,+oo) (2) 
= uo(x) Wx G (a, b) 

and its numerical resolution. In this paper we consider only the Euler for- 
ward in time - centered in space (EFC) schema (see jl]) for Finite Differences 
method, but for other numerical schemas one can apply the same consider- 
ations with analogous, but more complicated, demonstrations. Let h > a 
linear step for a partition (xj)j, < j < P, P = 2=2 ( of the interval [a, b] 
such that Xq — O, Xp — b, and let At the time step of the numerical schema. 
We denote by u n = (u™)j the numerical solution of problem (J5J) computed at 
time step n-th, hence at instant t = nAt. 

In EFC schema the values at time step n + 1 are so computed (see 

< +1 = < - - U U)> o< 3 <p-i (3) 

where F'(uj) is the value of j£(u) obtained by substitution of u with the 
discretized value u". We suppose that the schema (jHJ) satisfies the Courant- 
Friedrichs-Lewy condition (CFL) (see |3]): 

| F'(u]) | < 1 Vj G Z, Vn e N, V/i G M+ (4) 

This condition is usually admitted in the literature about numerical resolu- 
tion of PDEs and is related to the question of stability of solutions. 



Let iV e N, iV > 1, r e N+, k = ~, dt = ^r, u N the numerical solution of 
(|2*J) at time step with linear step h, u M the solution at time step M = Nr 
with linear step k. Note that CFL condition is satisfied in the second case 
too, because x = T- ^ ne som ti° n uM is computed on a grid finer than that 
of u N , being h > k, and hence it can be considered a better approximation of 
the exact solution, if exists, of the problem (0) . But u M has the disadvantage 
of a much greater computational cost, which can be very high for r >> 1 
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(see 0). 



We construct a piecewise polynomial interpolation v of u on the nodes 
of the partition with grid step h. Function v can be considered as a global 
continuous and differentiable representation of the discrete values u N on the 
interval [a, b]. Then one can compute, for every couple (xj,Xj+i) of nodes, 
the interpolating polynomial Vj = Vj(t) on a set of r values of t. In the 
following, we present an estimate of the error \vj(t) — u!f\, where t = —r 2 - 
(see Section 3), with the aim of establishing if v can be considered a good 
approximation for the solution of the problem (J2J). 

2 Construction of interpolating polynomials 

Let (xj,Xj+i) a couple of nodes internal to [a, b], hence 1 < j < P — 2. 
For simplicity we don't consider the boundary nodes xq and Xp, but with a 
suitable choice of a discretization formula for the spatial derivative |^ (see 
|Ej) w e could apply analogous considerations. We write, for more simple 
notation, p 1 = xj, p 2 = Xj+i, d\ = u^, d 2 = Uj+i- We construct a cubic 
polynomial q = q(t) = at 3 + bt 2 + ct + d, < t < 1, imposing the following 
four conditions: 



where q'(t) = 3at 2 + 2bt + c. Then, solving the following linear system for 
the unknowns a, b, c, d 



\ 3a + 2b + c = d 2 
one obtains the expression for the derivative of each cubic: 

q'(t) = (6p! - 6p 2 + 3di + 3d 2 )t 2 + (-6pi + 6p 2 - U x - 2d 2 )t + d x (6) 



1. q{0) 

2. o(l) 

3. g'(0) 

i g'(i) 



4; 




d = Pi 

a + b + c + d = p 2 
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3 Auxiliary technical results 



Let w = u N the numerical solution of problem (J2J) obtained with a EFC 
schema with grid step h (/i-grid) and initial condition it sampled using step h. 
We consider a couple of internal nodes (xj, Xj+i) on /i-grid, and a subpartition 
( x j,i)o<i<r, where r G N+, mod(r,2) = 0, = a^o, = a?^. Let 
f = = q'j(t), where qj is the cubic computed solving the system (JHJ). We 
call u an interpolant solution of the problem (J2J). We define the set t m — —, 
where m G N, < m < r. If s = X ° +1 ~ X: > ; we have x m = x Jim = ras + x^, 
therefore 

1 

t m = — = —T = ™~ -T = ~r = JT^- 7 

r r h ran h 

Let k = ~, and let u = u M the solution obtained with grid step k and 
uq sampled using k. For simplicity we denote u m the value of u on the node 
x m = ms + Xj, where j, as before, is index of the /i-grid. In the following we 
demonstrate some preliminary results for estimate the error \v m — u m \, where 

V m = V(t m ). 



Proposition 1 Ift G [0, 1], exists a value s G [0, 1] such that 

* = t~{ (8) 

Proof. If y = t 2 — t, then y < and y' — 2t — 1, therefore £ = \ is the so- 
lution of y'(t) = 0. At to the function y has a minimum, because y" = 2 > 0. 
The value is = — |. Therefore we have — ^ < £ 2 — £ < 0, from which (|8jl 
follows with s6 [0,1]. □ 

Proposition 2 For every m such that < m < r: 

3 3 
\v m - v r \ < -h + 2\di - d 2 \ + -\di + d 2 \ (9) 
2 4 

Proof. From © we have v r = v(l) = d 2 , therefore 



v m - iv | < |(6pi - 6p 2 + 3di + 3d 2 )£ 2 + (— 6pi + 6p 2 - 4rfi - 2d 2 )£ + d 1 - d 2 
< |(6p! - 6p 2 )(£ 2 - £)| + | (3^ + 3d 2 )£ 2 - (4d a + 2d 2 )£| + K - d 2 
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Then, using (JSJ), we obtain 

3 3 
\v m - v r \ < -\pi -p 2 | + 2\d x - d 2 \ + -\di + d 2 \ (10) 

and hence 0, because \pi — p 2 \ = h. □ 

In the same manner, using the fact that vq = v(0) = di, one proofs the 
following 



Proposition 3 For every m such that < m < r: 

3 3 
\v m - v \ < + 2\di - d 2 \ + -\di + d 2 \ (11) 

In the following propositions e is a fixed positive real number, w° is the 
initial condition uq of (j2J) sampled on the nodes of /i-grid and u° the initial 
condition sampled on /c-grid nodes. 



Proposition 4 // |u? - < ^ Vi, then \u? - uf +1 \ < Vz, 

Vn<M,nGN + . 

Proof. By induction. Let n = 1. We have 

H - u\ +1 \ < \u] - «?| + | u ° - + \ u ° i+1 - u\ +1 \ (12) 

Using the EFC schema (j3J), the CFL condition (jlj) and the hypothesis, we 
have 



<3 f 



3 A/ 3(M-1) 

that is the first inductive step. The intermediate inductive step is |w™ - 
<+il < ^w=zj- Then 

|< +1 - tijft 1 ! < |< +1 - <| + |< - < +1 | + \K+i ~ <ti\ 
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1, „ „ , , „ „ , 1 



— 2 1 M-l l—L I 1 

\\< +l - K\ + uu\ + 1< - < +1 | + \\ 



i+l 



<3- f 



1 + 


-l?/ n 
2 1 i+2 


-<l 








1 + 







3(M-n) 3 (A/-(n+l)) 

and this is the final inductive step. □ 

Proposition 5 // \v% - u° i+1 \ < pV Vi, i/ien - w™| < 3(JV /_ n) Vz, 

Vn < M. 

Proof. Let n > 0. Then 

l< +1 -<l< ~KVi-<-il 
<^(l< + i-<i + K-<-il) 

and hence, from Proposition (jlj) 

K +1 -<I<^T (13) 



3 (M-n) 



that is the thesis. □ 



In the next proposition, the index i is associated to the node Xi in the 
A;-grid, while the index j is associated to the node Xj in the /i-grid. 



Proposition 6 If \u® — < pr Vi, t/ien \i£ — Uj +1 \ < pr Vj. 

Proof. Using the fact that r < 3 r and M = Nr > (N + r) if r, N > 1, we 
have 

r— 1 i — 1 r— 1 

l I ^ ^ I I ^ ^ ^ ^ ^ ^ ^ 

|«j - < \ u j+f - u j+^ \ - - 1^ ^ = 3IV 

m=0 m=0 m=0 

therefore the proposition is verified. □ 

Note that by definition w® = vPa, hence from Proposition (jHJ follows that 
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K-<il<^ (14) 

Then the Propositions (JJJ) and (jSJ) hold for the solutions (w k ) too, using 
N instead of M: 



Corollary 1 If\u°-u° +l \ < pV Vi, tfien (w^-u^l < ^sy, |«#-tu* +1 | < 
^by Vj, V£;<iV. 

In particular, for k = N, writing Wj = d\ and Wj +1 = d 2 , follows that 

Corollary 2 If \u® — < pV Vi, t/ien |di — d 2 | < e- 

In the next proposition the simplified symbol uq = Uj 4 is the value of 
solution u at node Xj r of the &-grid, i.e. the node Xj of the /i-grid. 



Proposition 7 If \iq — u° +1 | < pV Vi, t/ien |v — u \ < e. 

Proof. By definition, vq = d\ = Wj = , where is the solution at 
time step N of (J2J) with initial condition u> . Then, using a telescopic sum 
like (a N - a ) = Y.k=i( a k ~ Ofc-i): 



AT M 

1^0 "Mo I 



fc=l h=l 

By definition, = hence the second term of the second member is null. 
For the first sum, from hypothesis and from Corollary (Q) follows 

N N 
k=l k=l 



N 



e 3^ - 1 e 

< 



3(at+i) ^ 3(Ar+i) 2 - 2 



For the second sum the computation, using the Proposition (jSj), is the same. 
Therefore the proposition is verified. □ 
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Using the simplified symbol u r = ufj +1 \ r for the value of u at the node 
xy+iv at time step M, in the same manner we can proof the following 

Proposition 8 If |u° — < pV Vi, t/ien |iv — w r | < e. 

In the next Proposition we remember that u m is the value of u = u on 
the node x m = ms + Xj. 

Proposition 9 Let 1 < m < |. — < p> Vz, t/ien |wo — w m | < me. 
Proof. Using Proposition (j^J) with n = M, we have 

m 
k=l 

Proposition 10 Let | < m < r. If \u® — i^+il < pY Vz, i/ien |w m — u r | < 
(r — m)e. 

Proof. Using Proposition (j^J) with n = M, we have 

r— m 

\u m - u r \ < \Um+k-i - u m+k\ <{r- m)e □ 

k=l 

4 The main results 

Now we can stated the main results on the estimate about the accuracy of 
the interpolant solution v. 

Theorem 1 Let e > 0, m G N+, 1 < m < r. If \u® — < pr Vz, then 

3 3 

\v m — u m \ < -h + -\di + d 2 \ + (min[m, r — m] + 3)e (16) 
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Proof. If m < |, let p = 0, else p = r. Using Propositions ((21) or (jSJ), (|Z|) 
or ©, © or {Till and Corollary (dJ), we have 

I | ^ | ^p I "I - I ^p ^p I I ^p | 

3 3 

< -h + 2e + — | g?i + c^2 1 + {min[m, r — m] + l)e 
because if m < | then mm[m, r — m] = m, else rranfm, r — m] = r — m. □ 

Corollary 3 Let e > 0, m G N +; 1 < m < r. There exists 5 > suc/i t/iat 
i/ |u° — < S Vi, then 

3 3 

\v m - u m \ < -h + -\di + d 2 \ + e (17) 
Proof. Let a = — — - — !—=•. Then u > 0, and from Theorem JT1), defined 

J ~ min[m,i — mJ+3 r ' *~~ » ' 

5 = pr, follows that if |w° — < 5, then 

3 3 

\v m — u m \ < —h + -\dx + d 2 \ + (mm[m, r — m] + 3)/i (18) 
But (min[m, r — m] + 3)/x = e, so the thesis is verified. □ 

The estimate stated in Corollary (JHJ) is of the same kind of the esti- 
mates established, in the case of finite element method, by Kuznetsov and 
Tadmor-Tang for approximated numerical solutions to scalar conservation 
laws equations (see P). 

In particular, inequality (fTT|) shows that the accuracy, respect to the fine u 
solution computed using a grid step k, of the interpolant solution v, com- 
puted using a grid step h = rk, is of the first order on h. In (|17|) the role 
of the time steps N and M is hidden in the third addendum e, as stated by 
Corollary ©. 

The next theorem gives a sufficient condition on the initial function uq for 
the validity of the hypothesis in Theorem (Q). 



Theorem 2 Let e > ? h > 0, N G N, N > 1. If u G C°([a,b]), then exists 
r G IR + such that, if < m < r, 

3 3 

\v m ~ u m \ < -h + -\di + d 2 \ + e (19) 
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Proof. The function u , continuous on a closed and limited subset of R, 
is uniformly continuous, therefore exists 8 > such that if x±,x 2 G [a, 6], 

— 2^2 1 < 5, then \uq(xi) — uq (x 2 ) \ < e. We can choose the parameter 
r G N+, and hence the time step M = rN, such that < 8. Therefore, 
with this choice of r, it can be applied the theorem ([TJ and, with the right 
choice of 8, the Corollary (j3J). □ 

In the next section we discuss the role of the term \d% + d 2 \ in some 
particular cases. 

5 Applications to particular cases 

As first application of the previous estimate, let F(u) = au, with a G R, 

a / 0. Hence we suppose that F is linear. In this case the numerical 

i 

schema EFC is stable in the norm \u\ 2 = yh^2j \ u j\ 2 ^j 2 until time step N, 

i.e. \u\ 2 < C n \uq\ 2 with C N = C(N) G R+, if At < (^) 2 ; the constant C N is 

equal to e 1 ^ (see jlj), while j is the index on /i-grid. If we choose the grid 
step h such that h < |a|, than the CFL condition is satisfied because 

. . A* h 

oh-<i-r<l 20 
h \a\ 

Lemma 1 Let m G N, m > 1, and (Cj)o<i<m Q> set of positive real numbers. 
If A, B G (Cj)j, then 

tm \ 2 

Proo/. We have A 2 < J2T=o C j and B * < E7=o C j, therefore 

(m \ 2 / m \ 2 

hence the thesis is verified. □ 

If uG C°([a,6]), let \u\oc = sup xe[ajb] \u(x)\. 
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Corollary 4 Let e > 0, N G N } N > I, h > 0, a / 0, h < \a\, At < (£) 2 
and uq G C°([a, b]). Then exists r G M+ such that, if < m < r, 



< 



\ u o\ 



+ e 



(23) 



Proof. From theorem (j2J) we have 

3 3 

\v m ~ u m \ < ~h + -\di + d 2 \ + e (24) 

where di = and d 2 = for a generic node Xj of /i-grid. Let s = ^ 

and | Wo |oo = su Pj\ w o,j\- From Lemma (0) and from stability in the norm | • I2 
follows that 



§14+41 <!(Ki + K.i)<§(E< 



0=0 




N O NAt 

-\w 2 < —re 2 Fob 



2/i2 2A 

3 I iVAt 

< - (s+ l) 2 e 2 I^qI^ 
and the thesis is verified because |u>o|oo < l^oloo- D 

The Corollary Q shows that the error \v m — u m \ is uniformly limited 
on interval [a, b] for every time step n < N. The presence of some kind of 
norm on Uq, like |«o|oo5 m the upper limitation is typical in literature about 
accuracy of numerical solutions (see jHj). 

As second application, let w the numerical solution at time step iV of the 
problem (J2J) on the /i-grid, and let Xj, Xj + % two nodes of this grid. We call 
the interval [xj, x,+i] an interval of rising turbulence at time step N for the 
/i-grid if (wjWj + i) < 0. The definition is based on the intuitive fact that a 
change of sign for the velocity field of a linear flow induces a superposition of 
particles in the motion of a fluid. In one dimensional case, in this situation 
the velocity u in the problem (j2J) might be not a classic function. 
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Corollary 5 Let e > 0, N e N, iV > 1, h > 0, u G C°([a,b}) and 
[xj , Xj + i] C (a, 6) an interval of rising turbulence at time step N for the 
h-grid. Let w N the numerical solution at time step N. Then exists r e M + 
such that, if0<m<r, v Q = w 1 ^ and v r = wf +1 , 

\v m -u m \ < -(3/i + 5e) (25) 
where u is the solution on k-grid, k = rh. 

Proof. From theorem (J2J we have 

3 3 

\v m ~ u m \ < -h + -\wf + wf +1 \ + e (26) 

From hypothesis we suppose Wj = < and Wj+i = w^ +1 > 0. Then exists 
w > such that Wj = —wwj+i. Theorem (j2J) holds, and from Corollary ((TJ) 
we have \vjj — u>j+i\ < e. Therefore 

\vjj — Wj + i\ = | — uiWj + i — uij + i \ = vjj + i\ — w — 1\ = Wj + i(l + w) < e (27) 
and hence tUj+i < < e. Also 

\wj — Wj+i\ = —Wj + Wj + i < e (28) 
hence \vjj\ = —Wj < e — Wj + i < e. Therefore we have 

\ w j + ^j'+il — \ w j \ + l^i+il — 2e (29) 
and from (j26|l follows that 

3 3 

\v m - u m \ < -h + -e + e (30) 

The thesis is verified. □ 

From the proof of the Corollary we may note that in an interval of rising 
turbolence the values of the solution at extremes are small. The particular 
limit case Wj = Wj+i = 0, for which, from theorem (J2J), \v m — u m \ < |/i + e, 
might be interesting for searching periodic solution u for which Uj = Uj+i = 
on the nodes of the /i-grid: 
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